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ABSTRACT 

In order to construct accurate point sources simulations at the frequencies relevant to 21 -cm 
experiments, the angular correlation of radio sources must be taken into account. This paper 
presents a measurement of angular two-point correlation function, w{6), at 232 MHz from 
the MIYUN survey - tentative measurements of w{9) are also performed at 151 MHz. It is 
found that double power law with shape w{d) = A9~'^ fits the 232 MHz data well. For the 
angular lenght of 0.2° ^ 6 < 0.6°, 7 « —1.12, and this value of slope is independent of 
the flux-density threshold; while for angular lenghts much greater than 0.6, 7 has a shallower 
value of about —0.16. By comparing the results of this paper with previous measurements of 
w{d), it is discussed how w{9) changes with the change of frequency and completness limit. 
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1 INTRODUCTION 

There has been considerable recent interest in using the highly red- 
shifted hyperfine line from Hi (21 -cm li ne) for astrophysical an d 
cosmological studies at redshifts z > 6 jPritchard & Loeb 2008h . 
At these redshifts, the 21-cm line is redshifted to metre wave- 
lengths, and extracting the cosmological signature requires ac- 
curate modeling and re moval of the fore g round Galactic and 
extragalactic em ission jSantos et al. 20051 [Morales et al. 200^ : 
IWang et al. 200^ . 

The diffuse Galac tic foreground emission fluct uates mainly 
on large angular scales l lde Oliveira-Costa et al. 2m§) . or on scales 
that are much larger than the expected angular fluctuations of the 
21-cm signal. Point source contamination from discrete extragalac- 
tic radio source, on the other hand, affects mainly small angular 
scales and potentially could be more problematic. A number of sur- 
veys of radio sources have been performe d at frequencies rel evant 
to the 21-cm tomography - see Figure 1 in dCohen et al. 2003h : and 
analysis of these catalogs have helped to bring some understanding 
about their statistical properties: for instance, it is known that the 
distribution of radio sources obeys a Poisson statistics with an ob- 
served angular clustering - see Table [H 

In this respect, our knowledge of the extragalactic radio source 
clustering properties is important because, if there is clustering 
on small angular scales, it could contribute power to a power 



The range of frequencies we consider in this analysis coirespond to a 
redshift range z Ri 6-8, which are well within the range being considered 
for the Epoch of Reionization. 




Figure 1. Footprints of the 6C 151 MHz catalogue (top, left), the 7C 
151 MHz catalogue (top, right), the MRT 151 MHz catalogue (bottom, left) 
and the MIYUN 232 MHz catalogue (bottom, right). All footprints were 
computed using the HEALPix projection I Gorski et al. 2005), on which lo- 
cations covered by the catalogs are in a brighter color (we adopted the con- 
vention in which an increase in flux density S corresponds to an increase in 
color's brightness). All catalogues are plotted in the interval of ^ 5 ^ 1 
Jy, and in Galactic coordinates with the Galactic center at the origin and 
longitude increasing to the left. 



spectral analysis; as a result, the clustering signal could be con- 
fused with the 21-cm signal that is being sought. Further, in or- 
der to construct accurate simulations at metre wavelength^ the 

^ Some aspects of both experimental design optimization and ac- 
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angular correlation of radio sources must be taken into account 
jOonzalez-Nuevo et al. 2003) . At the faint flux densities relevant 
for 21 -cm cosmological studies, the relative importance of the clus- 
tering contribution increases and could become an important con- 
tribution to power spectral analyses. 

This paper presents a measurement of angular two-point cor- 
relation function w{6) at 232 MHz and a tentative measurement 
at 151 MHz. By comparing the results of this paper with previ- 
ous measurements of 10(0), we assess how w{9) changes with the 
change of frequency and completness limit. In Section |2l we de- 
scribe the statistical tools as well as the surveys used in this analy- 
sis. The results and conclusions are presented in Sections [3] and [4] 
respectively. 



2 DATA ANALYSIS TOOLS 

2.1 The Angular 2-point Correlation Function 

The clustering of astronomical sources is quantified using the an- 
gular two-point correlation function w{9). One way to estimate this 
function is to compare the distribution of the objects in the real ca- 
talogue to the distribution of points in a random Poisson distributed 
catalogue with the same boundaries cHamilton 1993.) . or 



00(61) * RR(6') 



(1) 



where DD(6l), RR(6l) and DR(6l) are the numbers of data-data, 
random-random and data-random pairs separated by the distance 
6* + 561. 

The estimation of RR(S) and DR(6') requires a catalogue of 
objects distributed uniformly over an area with the same angular 
boundaries as the data catalogue. In order to produce such cata- 
logues, we used the "Sphere Point Picking Algorithm'lf] to gene- 
rate random cartesian vectors equally distributed on the surface 
of a unit sphere (to avoid having vectors "bunched" around the 
poles, as would happen if spherical co ordinates were used). See 
dde Oliveira-Costa & Capodilu po 200^ for more details on how 
these catalogues were produced. 

Incomplete knowledge of the completeness limit of a sur- 
vey can affect our determination of the extent of clustering 
jMagliocchetti et aL 1998). If a completeness limit for a given sur- 
vey cannot be found in the literature, we estimate a ballpark value 
from the survey's differential source counts dN/dS. We derive 
dN/dS by binning the sources in flux density with bins of 0.1 in 
width, and the bins are not weighted by S~^'^. 



2.2 The Catalogues 

All the catalogues decribed in this subsection are avaiable at 
\http ://vizier. cfa. harvard. edu/viz-bin/VizieR\ and are shown in Fi- 
gure [T] Below we provide a brief description of the data sets used 
in this analysis. 



tual data analysis require full-blown simulations of the sky sig- 
nal and knowledge about how it propagates through the instrument 
and the data analysis pipeline. End-to-end simulations are important 
for 21-cm experiments because of the many complicated issues re- 
lated to instrumental performance, ionospheric tu r bulence corrections , 
etc. /Morales & Hewitt 2004; B owman etal. 200d ; [Morales et al. 20051 ; 
iBowm an 2007; Bowman et al. 200^^ 
^ [http;//mathworld.wolfram.com/SpherePointPicking.html| 
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Figure 2. The measured source counts dN/ dS (top) and angular coiTela- 
tion function w{6) (middle, bottom) of the 6C 151 MHz catalogue. Top: 
The source count of the 6C-II region is plotted using the peak amplitudes. 
The red line is a single power law fit to the histogram (black line), for which 
dN/dS = 2.61 - 1.545. Middle: The angular correlation function of the 
6C-II region calculated for a flux density limit of 5 = 500 mjy. Bottom: 
The angular correlation function of the (6C-II + 6C-IV) region calculated 
for a flux density limit oi S = 500 mJy, restricting to Galactic latitudes 
\b\ > 5°. The yellow shaded regions in the middle and bottom panels re- 
present w{0) calculated using solely mock catalogues. 



The e**" Cambridge (6C) survey produced a catalogue of radio 
sources at 151 MHz with 34,418 discrete sources at an angular re- 
solution of 4.2' X 4.2' csc(5) . The data product is a set of seven sec- 
tions (named 6C-I, 6C-II, 6C-III, 6C-IV, 6C-Va, 6C-Vb, and 6C- 
VI) that maps the sky north of +30° declination, with limiting flux 
densities between 130 mJy and 200 mJy, depending on the section 
analvsed tealdwin et al. 1985l:lHa les et al. 1988l:lHales et al. 19"90l : 
iHales et al. 199 ll : iHales et al. 19 93a; Hale s et al. 1993bl) . The im- 
certainties on the positions and flux densities of many sources 
in the regions 6C-Va, 6C-Vb and 6C -VI are not well quantified 
Males et al. 1993iJ : I Hales et al. 1993bh . so we exclude these re- 
gions from our analysis. 

The Cambridge (7C) survey produced a catalogue of radio 
sources at 151 MHz with 43,683 discrete sources at an angular re- 
solution of 1.2' X 1.2' csc((5). This survey is composed by combin- 
ing 96 individual images at declinations greater than +21°. These 
images have completness limits between 120 mJy and 770 mJy, 
depending on the image analysed jHales et al. 2007h . 

The Mauritius Radio Telescope (MRT) is a Fourier synthe- 
sis array that has produced images of the sky covering the r egion 
IS*" < a < 24'' and -75° < S < -10° jNavak et al. 2009h . The 
resulting catalogue contains 2,784 discrete sources at an angular 
resolution of 4.6' x 4.6' with a completness limiQof about 1 ly. 

The MIYUN 232 MHz survey mapped the sky north of decli- 
nation +30° at an angular resolution of 3.8' x 3.8' csc((5), with 
an average noise level of 50 mJy. The principal data product 



Shankar 2009, private cominication. 
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7C. S„„ = 270mJy 




Figure 3. Measured angulai' con'elation function w{8) of the 7C (top) and 
MRT (bottom) 151 MHz catalogues. The angular con'elations are calculated 
at the flux density limit of 5 = 270 mJy (7C) and S = 1 Jy (MRT). The 
yellow shaded region is w{d) calculated using solely mocks. 



of this survey is a catalogue containing 34,426 discrete sourc es 
jZhang et al. 19971) . complete at the 250 mJy level jZhang 199^ . 







MIYUN, Sii^ = 250mJy. |b|>10° 
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3 RESULTS 

3.1 The catalogues at 151 MHz: 6C, 7C, & MRT 

The completeness limits of the the 6C-I, 6C-II, 6C-III, and 6C- 
rV regions were not pubhshed. Using the source counts dN/dS, 
as described in i]2.1l we estimate completeness limits of 200 mJy, 
500 mJy, 250 mJy, and 500 mJy, respectively. Figure |2] illustrates 
our estimation procedure. The top panel shows the differential 
source count (peak amplitudes) of the 8,275 objects in the 6C- 
II region. The red line corresponds to a single power law fit to 
the dN/dS distribution, for which dN/dS = 2.61 - 1.545. At 
flux densities below 800 mJy (vertical black dotted line), the num- 
ber counts begin to flatten, and the lack of faint objects becomes 
very important below 500 mJy (vertical black dashed line). We 
conclude that this survey region is incomplete at flux densities be- 
low 500 mJy. 

Figure |2] (middle) shows the measured w{0) of the 6C-I1 re- 
gion for our assumed completeness limit of 500 mJy. The mea- 
sured correlations are represented by the black squaresQi and dis- 
tances between data and/or random sources are measured in bins of 
0.09° . At this flux density limit, the 6C-II catalogue contains 4, 1 1 1 
objects. The results are consistent with zero. Similar results were 
obtained for the three remaining regions 6C-1, 6C-1II, and 6C-1V. 
We also investigated if w{9) changes with a change in bin size, in 
Galactic latitude cut, or flux density limit. There is no indication 
that modest changes in any of these quantities affect the results. 

^ In order to avoid having Galactic sources in our analysis, we adopted 
Galatic latitude limits and discarded sources at lower Galactic latitudes be- 
fore measuring w{9). 100 mock catalogues are constructed using the pro- 
cedure described in i|2. II with flux densities above the sensitivity limit of 
the data catalogue and a chosen Galactic latitude cut (if applicable). By 
cross-correlating the data with the 100 mocks, a set of normally distributed 
estimates of the congelation function is produced. The mean and the stan- 
dard deviation of this distribution are used as a value for the estimate and 
its uncertainty in the measurement of w{6) at each 8. The estimate (mean) 
and its uncertainty (the standard deviation) are shown as the black squares 
and their uncertainties (e.g., Figure|2] middle). Similarly, the 100 mocks are 
correlated with themselves. This result corresponds to the yellow shaded 
region shown in the correlation figures, and, as expected in a Poissonian 
distribution, w{6) is consistent with zero. 



Figure 4. The angular correlation function w{6) and the amplitudes A 
and power law indices 7 as a function of flux density determined from 
the 232 MHz MIYUN catalogue. Top: The angular correlation function is 
calculated for a flux density limit of 250 mJy, for a Galactic latitude cut 
of 10°. The red solid and the dashed blue lines are single power laws 
that fit the data, with w{9) = Ad Middle and Bottom: The ampli- 
tudes A and power law indices 7 are measured for various flux density 
limits from 250 mJy to 500 mJy. in steps of 50 mJy. The amplitude of clus- 
tering does not depend on flux density. 



As an extra test, w[9) was measured over grouped survey 
regions with similar completness limits. Figure |2] (bottom) shows 
w{6) calculated for the combined (6C-1I -1- 6C-1V) region, with a 
completeness limit of 500 mJy and a Galactic latitude cut of 5°. 
At this flux density limit, the (6C-1I -1- 6C-IV) catalogue contains 
6,806 objects, extending from -f 30° to +82° in declination. The 
data show the presence of a large scale correlation (which extends 
beyond 10°), but no signs of sma ll scale (9 < 1°) clustering. As 
discussed by j Blake & Wall 2003) . a possible explanation for this 
large-scale correlation is a varying areal density of sources on the 
sky, such as might occur from combining two 6C survey regions. 
In the case of 6C-II and 6C-IV, the former covers a declination 
range of +50° to +51° while the latter covers +67° to +82°. The 
change in projected interferometric baselines over this range of de- 
clinations would produce a varying surface brightness sensitivity, 
which in turn could affect the density of sources on the sky, and 
thereby spuriously enhance the measured value of 'w{9). 

Like the 6C catalogue, the 7C catalogue was extracted 
from a survey composed by multiple individual images that have 
completeness lim its varying between 120 mJy and 770 mJy 
jilaies et al. 2007h . Therefore,ui(6') was measured in each individ- 
ual region, as well as in aggrouped regions having the same com- 
pletness limit. Figure |3] (top) shows an example of the resulting 
angular correlation function for a region with a completeness limit 
of 270 mJy. At this flux density limit, the catalogue contains 2,872 
objects. Distances between data and/or random sources are mea- 
sured in bins of 0.09°. The results are consistent with zero. Similar 
results were obtained for the other regions in the 7C catalogue. We 
also investigated whether w{9) changed with a (modest) change in 
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Table 1. Published w(9)l values. 



Ref 


V 


A 


7 


w(e) 


Slim 




[GHz] 






[] 


[mJy] 


fde Oliveira-Costa & CaoodiluDO 2009) 


0.074 


0.103±0.026 


1.21 ±0.35 


0.2-0.6 


770 


This Work (6C) 


0.151 








500 


This Work (7C) 


0.151 








120-770 


This Work (MRT) 


0.151 








1000 


This Work (MIYUN, |b| > 10°) 


0.232 


0.096± 0.071 


1.12 ±0.11 


< 0.6 


250 


A X 10"=* 


fSeldner & Peebles 1981) 


0.178 






1.50-3.0 


3000 


(Renselink & Rotteerine 1999) 


0.325 


2.0 ± 0.5 


0.8 


< 1.0 


35 


(Blake et al. 2004) 


0.325 


1.01 ±0.35 


1.22 ±0.33 


> 0.2 


35 


(^bster 1977a) 


0.408 








250 


(Webster 1977b) 


0.408 








10 


fBlake et al. 2004) 


0.843 


2.04 ± 0.38 


1.24 ±0.16 


> 0.2 


10 


(Cress et al. 1996) 


1.400 


3.7 ±0.3 


1.06 ±0.03 


0.02-2.0 


3 


(Maaliocchetti et al. 1999) 


1.400 


2.68 ± 0.07 


1.52 ±0.06 


0.30-3.0 


3 


(Overzier et al. 2003) 


1.400 


1.2 ±0.1 


1.8 


> 0.3 


3 


fBlake et al. 2004) 


1.400 


1.49 ±0.15 


1.05 ±0.10 


> 0.3 


10 


fOverzier et al. 2003) 


1.400 


1.0 ± 0.2 


1.8 


> 0.3 


10 


fWebster 1977a) 


2.700 








350 


(Sicotte & Peebles 19^ 


4.850 






0.70-1.7 


45 


(Kpoiman et al. 1995) 


4.850 


4.01 


0.8 


0.30-1.9 


35 


fRenselink & Rottserine 1999) 


4.850 


6.5 ± 2.0 


0.8 


< 2.5 


35 


fLoan et al. 1997) 


4.850 


10.0 ± 5.0 


0.8 


< 2.0 


50 



^ w(6) is fitted by a power-law of the form Ad T 
Slim is the limiting flux density. 



bin size or Galactic latitude cut, and there is no indication that any 
of these changes affect the results. 

Figure |3] (bottom) shows the measurement of wiff) for the 
MRT catalogue, with a flux density limit of 1 Jy (black squares). 
At this flux limit, the catalogue contains 2,294 objects. Distances 
between data and/or random sources are measured in bins of 0.09° . 
The results are consistent with zero. We also investigated whether 
vj[Q) changed with a (modest) change in bin size, flux density limit, 
or Galactic latitude cut, and there is no indication that any of these 
changes affect the results. 



there is a small angular scale signal (maybe caused by clustering) 
added to a large angular scale signal in the surve)[f]. 

We also investigated whether w{Q) changed with a (modest) 
change in bin size or Galactic latitude cut. Variations in the bin size 
did not affect the results, but for Galactic latitude cuts greater than 
+ 15°, is consistent with zero. 

The angular correlation function w{Q) was also calculated for 
various flux density limits from 300 to 800 mJy, in increments 
of 50 mJy. Above the flux density limit of 500 mJy, the correlation 
for small angular scales (6* < 0.6°) approaches zero. As shown in 
Figure |4] (midd/e and bottom), for a 10° Galactic latitude cut, the 
amplitude of clustering does not depend on flux density. This same 
result was observed in previous angular correlation analy sis, e.g., 
telake et al. 2004lde Oliveira-Costa & Capodilupo 2009() . 



3.2 The MIYUN catalogue 

Figure|4](top) shows the determination of w(ff) for the flux density 
limit of 250 mJy (black squares). Distances between data and/or 
random sources are measured in bins of 0.09°. It is apparent that 
the angular correlation function cannot be fit with a single power 
law, but req uires two power laws, each of the form wiff) — A6^"' 
jPeebles 198Ci) . where A is a measure of the amplitude of the 
average enhancement of the number of radio sources at a parti- 
cular point in the sky. Fitting the data with a double power law 
model yields a weak correlation with A = 0.096 ± 0.071 and 
7 = -1.12 ± 0.11 (red sohd line, with x^=0.01), and A = 
0.236 ± 0.092 and 7 = -0.16 ± 0.05 (blue dotted line, with 
X^=0.24). The break in the angular correlation function occurs be- 
tween 0.4° ^ 6^ < 0.6°. This break could be the indication that 



4 CONCLUSIONS 

We have used existing 151 MHz and 232 MHz catalogues to 
estimate the angular correlation function 'w{9) relevant for both 
simulation and analysis of 21 -cm cosmological observations. 
At 151 MHz (corresponding to z ~ 8.4), our estimate of 'w{6) 
is determined from the 6C, 7C, and MRT surveys. In all cases, the 
results are consistent with zero, implying no observed clustering 
of radio sources, at flux density limits ranging from about 0.3 Jy 
to 1 Jy. At 232 MHz (corresponding to z ~ 5), we found that 

° As pointed out by < Blake & Wall 20021) . in a wide-area survey under- 
taken with an interferometer (such as the MIYUN), large scale gradients in 
completeness limit typically appear as a function of declination. 
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10(6') can be fit by a two broken power laws, each with shape 
w{6) — A6~'' , with a break at an angular scale of 6 ^ 0.5°. 
At small angular scales, with a Galactic latitude cut of 10° and a 
flux density limit of 250 mjy, we find A = 0.096 ± 0.071 and 
7 = —1.12 ± 0.11 (with x^=0.01), while, at large angular scales, 
we find A = 0.236 ±0.092 and 7 = -0.16 ±0.05 (with x^=0.24). 
The value of 7 at small angular scales is consistent with that mea- 
sured from radio catalogues at higher frequencies (TableO. 

From Table[T] it appears that the power law index 7 of the an- 
gular correlation function is essentially constant over the frequency 
range from 74 MHz to 232 MHz (20 < z < 5), and it is also 
consistent with that determined from higher frequencies, at least as 
high as 1400 MHz. Strikingly, however, it is not clear that the an- 
gular scales for which this power law index applies are consistent 
across the range of frequencies: at the lower frequencies, cluste- 
ring seems to be important on angular scales of 9 < 0.6°, while 
at higher frequencies, the clustering is important on angular scales 
6 > 0.3° . These are overlapping ranges, but not entirely consistent. 

It is importatnt to point out that the classes of sources being 
probed at the different frequencies are also not entirely the same. 
Source counts at 1400 MHz indicate that strong radio sources (e.g.., 
FR II radio galaxies) dominate at flux densities above about 10 mJy, 
while star forming galaxies become important at lower frequen- 
cies. The current flux density limits for surveys at frequencies 
around 200 MHz (and below) are about 500 mJy. Assuming a no- 
minal spectral index of a = —0.7 (S^ oc v"), these flux den- 
sity limits imply flux densities of about 125 mJy at 1400 MHz, 
which is well above the flux density limits probed by the surveys 
around 1000 MHz, and consistent with the notion that the current 
generation of low-frequency surveys are dominated by powerful 
radio sources. Consequently, one potential explanation is that the 
possible inconsistencies are due, in part, to the different source pop- 
ulations being probed. However, it is also notable that the clustering 
amplitude A is fairly marginal (always less than 4a, and in some 
cases consistent with zero) at frequencies lower than 1400 MHz; 
while at the higher frequencies, by contrast, the clustering ampli- 
tudes can exceed a significance of lOcr. 

For the purposes of 2I-cm cosmological observations, either 
simulating sky models or analyzing low radio frequency observa- 
tions, we conclude that it is acceptable to scale the angular clus- 
tering results from higher to lower frequencies. One caveat to this 
conclusion is if a strongly clustered, steep-spectrum population of 
objects exists. Using a fiducial 1400 MHz flux density of 10 mJy 
allows us to estimate how deep future low-frequency surveys might 
need to be in order to assess this possible inconsistency (or the ex- 
istence of a possible steep-spectrum population) in the scale of an- 
gular clustering. With the nominal spectral index of —0.7, we esti- 
mate that a survey around 150 MH4J would need to reach a limit- 
ing flux density of about 50 mJy (implying a thermal noise limit of 
about 7 mJy). 
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Such as the proposed Million Source Shallow Survey (MSSS) 
with the Low Frequency Array (LOFAR) - more information at 
|http://www! as tron.nl/general/lofar/lofar I 
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